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Abstract 

We use a four-dimensional lattice calculation of the full-QCD (quantum chromodynamics, the non-abliean gauge 
theory of the strong interactions of quarks and gluons) path integrals needed to determine the masses of the charmed 
and bottom baryons. In the charm sector, our results are in good agreement with experiment within our systematics, 
except for the spin-1/2 E ce , for which we found the isospin-averaged mass to be H ee to be 3665 ± 17 + 14*2„ MeV. 
We predict the mass of the (isospin-averaged) spin-1/2 Q. cc to be 3763 + 19 + 26+™ MeV. In the bottom sector, our 
results are also in agreement with experimental observations and other lattice calculations within our statistical and 
systematic errors. In particular, we find the mass of the Q.t to be consistent with the recent CDF measurement. We 
also predict the mass for the as yet unobserved "E' h to be 5955(27) MeV. 
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Experimental and theoretical studies of charmed and 
bottom hadrons have been the focus of vigorous re- 
search over the last several years. In particular, singly 
and doubly heavy baryon spectroscopy has received 
significant attention, mainly due to the recent exper- 
imental discoveries of both new charmed (SELEX) 
and bottom baryons by DO and CDF. In addition to 
these discoveries, there are still many states of heavy 
and doubly heavy baryons remaining to be discovered. 
The new Beijing Spectrometer (BES-III), a detector at 
the recently upgraded Beijing Electron Positron Col- 
lider (BEPCII), has great potential for accumulating 
large numbers of events to help us understand more 
about charmed hadrons. The antiProton ANnihilation at 
DArmstadt (PANDA) experiment, a GSI future project, 
and the LHCb are also expected to provide new results 
to help experimentally map out the heavy-baryon sec- 
tor. For these reasons, lattice quantum chromodynamics 
(QCD) calculations of the spectrum of heavy baryons 
are now very timely and will play a significant role in 
providing theoretical first-principles input to the exper- 
imental program. 

Lattice QCD is a discretized version of four- 
dimensional continuum QCD. It specializes in studying 
the strong-coupling regime of QCD, where perturbative 



approaches converge poorly. As in continuum QCD, we 
calculate an observable of interest through a path inte- 
gral: 
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where X.qcd is the sum of the pure-gauge and fermion 
Lagrangian, O is the operator that gives the correct 
quantum numbers for our observable, and Z is the par- 
tition function of the space-time integral of the QCD 
Lagrangian. It is straightforward to carry out this path 
integral numerically within a finite space-time volume 
and under an ultraviolet cutoff (the lattice spacing); in 
some physical cases, these approximations to the con- 
tinuum can dominate the systematic error; we must to 
carry out calculations using multiple volumes and lat- 
tice spacings to extrapolate to the continuum world. 

In this work we use the "coarse" lattice gauge 
configurations (with lattice spacing a as 0.125 fm) 
generated by the MILC Collaboration, in which the 
asqtad-improved Kogut-Susskind action is used for the 
fermions in the vacuum. However, for the valence light 
quarks (up, down and strange) we use five-dimensional 
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Shamir domain-wall fermion propagators with a fifth 
dimension of extent L5 = 16 and a mass parameter 
Ms = 1.7. We use hypercubic-smeared gauge links to 
minimize the residual chiral symmetry breaking. The 
pion masses range from 290 to 597 MeV, and the lattice 
volume is 20 3 x 64. 

For the charm quark we use Fermilab action, which 
controls discretization errors of 0((amg) n ). The inter- 
polating operators we use for the J — 1/2 singly and 
doubly charmed baryons are 

Ac : jfiifiCysqboZ, (h c : Cy 5 Q^q k u , 
Sc : e'HqilCysqbQi O n < : ^iqfCysQ^q], 
On x : e^Cfif Cy 5 «&& On a Cy 5 q{)Q k c , 

K : -^e* [tefCy,^)?* + (<?f Cy,^] , (2) 

where q u ,d,s are the up-, down- and strange-quark fields, 
and Q c is the charm-quark field. Given these actions 
and operators, we construct two-point functions using 
the path integral 

C k (t, t ) = J] {Ohix, t)Ol(x, t )) , (3) 

.? 

where Oh is an interpolating operator of the hadron h 
and (...) indicates an average over the lattice gauge en- 
semble. The correlation functions are calculated with 
gauge-invariant Gaussian-smeared sources and point 
sinks. The two-point correlators of the hadrons behave 
at large time like 

C h if) = Ae- E °' 7 (4) 

and the ground-state mass can be obtained from a 
single-exponential fit. 

We use static-quark action (with HYP-smeared links) 
for the bottom quarks. In the static limit, the bottom- 
quark propagator is reduced to a Wilson line. We can 
construct the static-light hadron correlation functions by 
contracting the static-quark line with light-quark prop- 
agators to create gauge-invariant correlation functions; 
the bottom-light baryon correlation functions have the 
form 

Grix, t; xq, t ) = (q" f ix, f)Tq b f ix, f)e cab 

P cc \x, f; xo, to)q~f,ixo, to)Tq b fi^o, foKw) , (5) 

where / e {u,d, s} is the light-quark flavor index and 
P cc is the Wilson line connecting the source and the 
sink which are separated by time t. Since the gamma 



matrix insertions in the baryon correlators are sym- 
metric iCy^) or anti-symmetric (Cys), the correspond- 
ing flavor wavefunctions are symmetric (sextet) or anti- 
symmetric (antitriplet), respectively. This flavor sym- 
metry structure is automatically enforced without ex- 
plicit symmetrization or antisymmetrization of the fla- 
vor wavefunctions. In the static limit the hyperfine 
splittings are exactly zero, so we cannot compute such 
spin splittings. For the physical bottom quark the hy- 
perfine splittings are less than one percent. When 
comparing to experiment, we identify our states with 
the lowest-spin physical particles. In the case of bot- 
tom hadrons, we study the mass splitting between two 
hadrons CA,(f)/C/n(0; this eliminates the unknown ad- 
ditive shift of the mass caused by using static (infinitely 
massive) quarks. 

Our results after extrapolation to the physical pion 
mass are summarized in Fig.Q] We also compare our re- 
sults with the dynamical lattice charmed-baryon calcu- 
lation j5t] and quenched calculations iSHHl, shown in 
Fig.Q] Ref. |H performed a calculation with staggered 
light quarks and an interpretation of Fermilab action 
defining the charm mass by the kinetic mass instead of 
the rest mass. They observe a significant lattice-spacing 
dependence in their calculation, which is very different 
from previous quenched calculations [7]. Their calcula- 
tion covers ensembles at more than one lattice spacing, 
but their results are not as precisely determined statisti- 
cally. Their uncertainties given here are merely the sta- 
tistical ones; no systematic uncertainties are given yet. 

We compare the doubly charmed baryons with pre- 
dictions from theoretical models; although the SELEX 
Collaboration has reported a first observation of dou- 
bly charmed baryons, searches by the BaBar, Belle and 
Focus Collaborations have not confirmed their results. 
This makes it interesting to look back to theory to see 
where the various predictions lie. We compute the mass 
of E cc to be 3665 + 17 ± 14+° g MeV, which is higher 
than what SELEX observed. 

We compare with a selection of other theoretical re- 
sults, such as a recent quark-model calculation, rela- 
tivistic three-quark model, the relativistic quark model, 
the heavy quark effective theory and the Feynman- 
Hellmann theorem. Most theoretical results suggest that 
the S cc is about 100-200 MeV higher. The Q. cc mass 
prediction made by this work is 3763 + 19 ±26 *lg MeV, 
and the overall theoretical expectation is for the Q. cc to 
be 3650-3850 MeV. We hope that upcoming experi- 
ments will be able to resolve these mysteries of dou- 
bly charmed baryons. In the future we plan to extend 
these calculations to a second lattice spacing, address 
systematics due to continuum extrapolation, and extend 
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Figure 1: (Left) A summary of charmed baryon masses (in GeV) calculated using LQCD. Our main results [1] are the lighter (orange) points 
taken from splitting extrapolation and the darker (blue) points taken from direct mass extrapolation. (Right) Comparison of mass splittings with 
all available 2+1-flavor lattice calculations of bottom baryons. The square (blue) points are extrapolated using the fl-mass reference scale Q]; the 
solid error bars indicate the statistical error, and the dashed bars indicate the total errors (including the estimated systematics). The solid (red) bars 
indicate the experimental values given in the PDG, where available. For the Cl b , we show both the DO result 1 3] (upper-right, purple) and the CDF 
result |4] (lower-left, magenta). 



our spectroscopy to spin-3/2 states. 

Our results for bottom baryons are shown on right- 
hand-side of Fig. [T| along with other 2+1 -flavor 
works 

All four calculations include the split- 
tings Efc - A h and E/, - A/,; we calculate the additional 
splittings E' b - £/, and Q/, - E' h in order to make direct 
comparisons. We again see good agreement amongst all 
lattice calculations, and mild scaling in the Na et al. re- 
sults. Agreement with experimental values as given by 
the PDG is fairly good. The discrepancy in the £/, - A;, 
splitting may be due to discretization effects, as sug- 
gested by the trend visible in the Na et al. results. 

Our results agree well with experiment in the cases of 
known mass splittings. In the case of Q/, our calculation 
is in agreement with the recent CDF result and several 
standard deviations away from the DO result. Such a 
conclusion holds for all other lattice calculations of the 
Q. b mass. Our results for the E' h mass splittings are a 
prediction since a' b has not yet been observed. Using 
the splitting from the A;, and summing with the experi- 
mental value of that baryon, we predict the mass of the 
a' b to be 5955(27) MeV, where the error given combines 
in quadrature all our errors. For more details about the 
lattice calculation (including the chiral extrapolations 
to physical pion mass and systematic error estimations) 
and references, we refer readers to Refs. 

In the future we plan to extend these calculations us- 
ing the anisotropic gauge ensembles generated by the 
Hadron Spectrum Collaboration (HSOlU d. The 
anisotropy (giving finer temporal lattice spacing while 
keeping spatial one large enough to keep the finite- 
volume systematic error small) will provide finer res- 
olution not only to improve the signal quality but also 
to reliably probe the excited states. In the charm sector, 
we will extend the calculations to J — 3/2 baryons and 



adopt a relativistic improved action for the bottom quark 
to get rid of the mass redundancy in the bottom sector. 
A recent paper by HSC 11311 demonstrates a new tech- 
nique, optimized for complicated operators, which may 
further improve the calculations and in particular allow 
for extraction of the excited-state spectrum of both the 
charm- and bottom-baryon sectors. 
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